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Abstract 

It is proved that true Goldstone bosons develop coherent fields whenever the 
associated charges of the matter particles are not conserved in a macroscopic 
scale. The sources of the Goldstone fields are the time rates of quantum 
number violation. The case of neutrino flavour oscillations is studied with 
application to Supernovae. It is shown that if the Lepton numbers break 
at the Fermi scale, the v potentials and oscillation patterns change in the 
periods of largest v fluxes. In this way, v e <-> oscillations may occur in 
the first instants of v emission. 

PACS numbers: 14.80.Mz, 14.60.Pq 
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I. INTRODUCTION 



In theories with spontaneous breaking of global symmetries the Goldstone theorem 
asserts the existence of massless particles - the Goldstone bosons - one per each non-zero 
quantum number of the vacuum. In principle such massless scalar bosons could mediate 
long range forces and lead to observable coherence effects just like the gravitational force. 
However, the very structure imposed by exact global symmetries at the Lagrangian level 
implies that these goldstons have nothing but derivative couplings with matter fields. It 
was shown that if the symmetry constraints are relaxed either through anomalies or with 
soft symmetry breaking terms then, scalar couplings of pseudo-goldstons with fermions 
naturally arise. Consequently, macroscopic forces are developed by matter aggregates but 
one price is paid: the intermediate bosons get non-zero masses which means that the 
range of those forces is finite. In this paper it is shown that under special conditions a 
true Nambu- Goldstone boson may develop a macroscopic field, the source being the rate 
of non-conservation of the matter quantum number associated with that field. That may 
be significant in astrophysical compact objects and the case of supernovae is analysed. 

To be definite I consider an abelian global symmetry U(1)a where A denotes gener- 
ically any linear combination of flavour lepton numbers such as L e — L M or the total 
lepton number L. After symmetry breaking non-invariant terms are generated in masses 
and/or couplings but the symmetry is still realized non-linearly and strongly constrains 
the couplings to the Goldstone boson 0. Under U(1)a, the matter fields and transform 
as 

ip a -> e - iaAa ip a , (la) 



-> + a V A , (lb) 

respectively, where Va is the scale of symmetry breaking. 

Consider for instance the lepton masses. Let z/£ denote all left-handed neutrino fields 
including the standard weak interacting partners v\,v^, v T L and any extra neutral singlets. 
The mass terms generated either radiatively or at tree-level are necessarily accompanied 
by couplings to completely determined by the lepton charges A a : 

C m = -m ab vf C v\ e *(A.+A»)*/VA + h.c. - m e £ £ , (2) 

where a sum over a, b and i — e, fi,r is assumed. Moreover, in any other term of the 
Lagrangian each field ip a always comes with a factor exp(z A a 0/VA). Consequently, a 
spatially uniform would be absorbed in the mass eigenstate fields leaving no couplings 
with 0. Matter only cares about the gradient of and that, as shown in |L], is better 
expressed in what I term as weak-unitary basis: 

i) a (old) = e - * Aa4,/VA i> a (uni) . (3) 

The new fields preserve the original quantum number assignments but are invariant under 
the global symmetry transformations ([!]). They reduce now to a translational invariance 
— ► + (3 and therefore, the goldston interactions can only depend on the derivatives. 



2 



In the unitary basis the effective Lagrangian writes as 



£ = £ mdi-Y.vf Cm ah v h L + h.c. 



e,Af,T a,b 



+ ld^d»<f> + J£d^/V A +--- 



(4) 



All of the goldston interactions are condensated in the current 



(5) 



which includes the matter current derived from the kinetic Lagrangian of ip (old) but also 
in general model dependent higher order corrections. The <fi equation of motion reads 



and clearly the source term is not some sort of matter density. 

The divergences of x/j^^i/j and ipj^j^ip only yield pseudo-scalar ipj^tp charges that 
do not generate long range T/r' fields, but rather spin dependent '1/r 3 ' potentials @. 
In the intent of escaping the derivative couplings law, models were made with explicit 
symmetry breaking terms. These may be anomaly terms [Q, as in the axion case 0, or 
soft symmetry breaking terms [BJ. In this way scalar couplings with the goldston field 
are produced but at the cost of getting a mass for the so-called pseudo Nambu-Goldstone 
boson which ultimately means a finite range for the long range forces. 

Here, I will explore the fact that the source term is different from zero whenever the A 
number of the matter particles is not conserved. To be accurate, there is a Noether current 
associated to the symmetry — > <p + a namely + V^d^cf), that is exactly conserved, but 
the non-standard interactions do not conserve the fermions A number as a result of the 
very spontaneous symmetry breaking. 

Non-conservation of may occur in macroscopic proportions in a stationary basis in 
stars, or in Supernovae for a short period of time. For instance, the total lepton number 
is violated at a constant rate in the process of neutrino emission of a star if u e has a 
Majorana mass: the non-conservation takes place in the neutrino production reaction 
where a fraction of the neutrino wave function is an anti- neutrino v e . This fraction is of 
course suppressed by the ratio m 2 u /El. That is not the case if a significant portion of 
the v e spectrum is converted into anti-neutrinos under appropriate magnetic field and 
transition magnetic moments. Another example is a Goldstone boson associated with the 
breaking of L e — L^. If there is resonant conversion of z/ e s into u^s then the goldston 
source has macroscopic dimensions and a long-range scalar field is generated. 



The solution of the <\> equation of motion is found using standard Green function 
methods. In the source term, d^J 1 ^ = pf represents the rate of creation of A-charge per 
unity of time and volume. In a stationary regime where this rate is constant or varies 
slowly enough one gets a Coulombian field 



(6) 



II. GOLDSTONE FIELDS 
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47T J Va \r — x\ 

where the bracket means expectation value of the operator over the appropriate particle 
state density. The Lagrangian of Eq. (f|) tells however that the fermions couple to the 
gradient = d^fi/V^. In the stationary limit A\ is negligible and the vector A^ = 
— V0/Va. obeys a Gauss law: 



V-A A = -(pa + V-J A )/V? 



It implies that the flux of Aa over a closed surface is proportional to the time rate Qa °f 
creation of A-charge inside the surface. In turn Qa * s equal to the difference between the 
flux of ( Ja) an d the flux of the current J° that would exist if there was no A violation 
inside the surface. Hence, 



A A = - {(Ja) - J?) K ■ (9) 

The value of J° is determined for a given p\ and obeys the identity p\ + V ■ Jl = 0. One 
reaches the same solution by noting that the Noether current + V^d^cj) is conserved. 

A typical situation of potential interest is the flavour violation driven by neutrino 
oscillations in their way out of a star. The quantum numbers are the partial lepton 
numbers Lg = L e , L^, L T or linear combinations of them and the bare currents are 

X = l-f t + ^ A + E WW* »l , (io) 

where z/£, a > r, denote any non-standard neutrinos. To be specific suppose that L T is 
spontaneously broken. The Goldstone boson <p T associated with it couples to a current 
Ja that is essentially j£ but contains in general higher order terms arising from radiative 
corrections. Their existence implies that even if L T is conserved in a certain process, the 
very reaction rate i.e., of particle production or absorption, is enough to yield a goldston 
field. However those terms are of higher order and I will limit myself to the bare tree-level 
currents. 

For definiteness I assume that all non-standard neutrinos are extra-heavy singlet neu- 
trinos that cannot be produced at energies of interest and in addition, no significant 
neutrino <-> anti-neutrino transitions occur either because the light neutrino masses are 
too small or because transition magnetic moments do not come into play. 

Suppose that in the way out of a star u e s convert into z/ T s and the flux of converted 
neutrinos is j(u e — ► u T ). Then, a macroscopic goldston field <p T exists with non-zero 
gradient outside the sphere where the conversion takes place. In the Supernovae case one 
also gets a flux j(u T — > u e ). Eq. gives A T = —5j T /V? with 

& = !{v e -> v T ) - ]{v T -> v e ) . (11) 

If the conversion involves anti-neutrinos one has instead Sj T = —j{p e v T ) + j{p T — > u e ), 
because the antineutrinos have negative Lepton numbers. 

The field <p T yields a potential energy for the r neutrinos. Generalising now to any 
quantum number Li and goldston field the interaction Lagrang ian jC< m t — Ji^Ag in 
Eq. (Hf) gives for the standard neutrinos V£, 
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(12) 



The potential energy of a neutrino with velocity v v is therefore 

V Vl = -A° e +A e -v u = -V Pe . (13) 

In the stationary limit A® = and for radially moving neutrinos V is just the radial com- 
ponent A\. The interesting point is that such kind of potential is not universal in flavour 
and in addition the sign of the potential is not solely determined by the matter content 
of the medium where neutrinos propagate as it also depends on the flavour structure of 
the neutrino fluxes. This allows in principle other than the standard MSW type of oscil- 
lations. The most interesting effect is probably the occurrence of u e <-> u x in a situation 
where it would not normally occur i.e., with a mass hierarchy m Ue < m v , m„ T . In the 
following I model such kind of phenomenon with application to Supernovae neutrinos. 

If a partial Lepton number is spontaneously broken it is natural that the same happens 
for all the three flavours. Let <7j be the set of scalar fields, singlets of SU(2)xU(l), with 
quantum numbers Aj = (L e , L M , L T \ under the symmetry group U(l) e x U(l) x U(l) r . I 
assume that the vacuum expectation values (<Tj) break all the Lepton numbers producing 3 
Goldstone bosons. By expanding around the vacuum as cr, = exp{— % £ A Aj}((<7j) + •••), 

A 

where the dots stand for the massive degrees of freedom, one obtains for the Lagrangian 
of the Goldstone fields 

£ = ^M^ A ^ M + iP/, (14) 

a sum over the quantum numbers A, M is implicit. The matrix 

VZ M = 2Y,A l M l \(cr t )\ 2 (15) 

is in general non-diagonal, the sole condition being that some of the fields <7j have simul- 
taneously more than one quantum number different from zero. Such mixing implies that 
the Goldstone bosons can mediate an action at a distance between two distinct currents 
as the equations of motion show: 

W A = -E G amW, (16) 

M 

where G\m is the inverse matrix of V£ M . Hence, in a stationary regime the vector field 
solutions are now 

Aa = -J2 G am((Jm)-4i) (17) 

M 

instead of Eq. (|9|). The expression of the one-particle potential remains the same as in Eq. 
(PI) because the interaction Lagrangian is the same as before namely, A^^ = j£ <9 M £ A . 

It is not the aim of this paper to construct specific models but as a proof of existence 
I present here a suitable adaptation of the abelian familon and seesaw Majoron 
models. Three neutral singlets N^N^N^ join the standard lepton iso-doublets L e ,L M , 
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L T , the charged singlets e^,/x^,r fi and the Higgs doublet $. The fields only have 
Yukawa interactions namely, flavour conserving Dirac couplings 

C D = -(m e L e $e R + tJ, e L e $N e R )/($ ) +h.c. , (18) 

plus Majorana couplings with a set of scalar singlets ae m = a m t- 

L m = -\H (K T CN%) M £m a em /(a em ) . (19) 

£,m 

The invariance under U(l) e x U(l) x U(1) T determines the (L e , L M , L T ) quantum numbers 
of the scalar fields namely, a ee (— 2, 0, 0), cr ejU (— 1, —1,0), 0^(0, —1, —1), etc. The expecta- 
tion values of these fields break all the three Lepton numbers which gives three Goldstone 
bosons. The matrix Ge m that enters in the goldston equations of motion, Eq. fll6D , has 
full generality with respect to its flavour structure. For Majorana masses M much higher 
than the Dirac masses fie, the neutrinos separate in heavy N e with mass matrix M and 
very small mixing with the three light neutrinos, ~ that acquire a Majorana mass 
matrix m tm = —fi^ji m Mi m . The magnitudes relevant for the subsequent application lie 
around m ~ 10 _2 eV, and M = 0.1 — 1 TeV, \i ~ 30KeV, very far from the original |7j 
seesaw GUT scales but not excluded from first principles. 

It is known that this type of singlet Majoron or familons easily escape present as- 
trophysical bounds on the couplings of Nambu- Goldstone bosons. Here, they couple 
primarily to neutrinos with pseudo-scalar couplings g ~ vn v jV^ which, for the specified 
range of neutrino masses, are far below 10~ 4 to play a role in Supernova collapse dynamics 
||, and even below the limit of ~ 10~ 8 ' 5 to produce Supernova cooling through singlet 
Majoron emission ||. On the other hand, pseudo-scalar couplings to electrons, that could 



be responsible for energy loss in stars ||10|| , only arise through radiative corrections and 
are so further suppressed. 



III. SUPERNOVA NEUTRINOS 

If the cause of the well known solar neutrino anomalies lies in neutrino flavour oscil- 
lations then one expects that the same phenomenum occurs with Supernovae neutrinos. 
In the accepted understanding of Supernova type-II all neutrino species including their 
anti-particles are emitted with similar importance after a first electron neutrino burst. As 
a working frame I adopt here the common prejudice that the neutrino mass eigenstates 
are in first approximation - but not exactly - V\ ~ v e , z/ 2 ~ u^, 1/3 ~ v T , with the following 
mass hierarchy: rri\ < m 2 < m 3 . In their way out from the star the electron neutrinos have 
a potential energy that exceeds by Vw the equal potentials of and v T . Vw = V2Gpn e 
[|TT| . p~2[1 is positive and proportional to the electron density n e and therefore the neutrinos 
propagate from the region of emission where Vw is much larger than m 3 (so I assume) 
to the outer space where Vw goes to zero. It is then in principle possible to have MSW 
resonant conversions from v e into v T and vice- versa |131J14|1. The anti-particle v e on the 
contrary, receives a negative potential — Vw and cannot be resonantly converted into the 
other flavours, given the assumed mass hierarchy. 
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In a Supernova the u e flux is larger than the v T flux and consequently more u e transform 
into u T than the reverse so making a net L e — L T violation. As a result, the goldston fields 
4>ei (j>(i, (fir develop as described above. Applying Eq. ( |iT| ) with flux transfers 



$Jt = S]e = j(Ve ~> Mr) ~ j( U T ~> f e ) , (20) 

one derives the vector potentials 

A e = (G ee -G eT )5j T , e = e,fi,T. (21) 

They vanish of course inside the sphere where v e <-> z/ r takes place. The magnitude 
of the constants G\m is directly given by the inverse square of the Lepton symmetry 
breaking scale, but the flavour structure is quite arbitrary and is not directly related to 
the neutrino masses as they also depend on the Yukawa couplings. The matrix G\m is 
only constrained to be symmetric and have positive eigenvalues as demands the goldston 
kinetic Lagrangian. 

An interesting case is the following hierarchy: \G ee \, |G eM |, \G eT \ « \G^ T \,G TT . Then, 
the vector A e is much smaller in magnitude than and A T . Far enough from the 
neutrinosphere and assuming spherical symmetry, the v potentials are well approximated 
by 

V Ve = V w + A e , \A e \«\V„J, 

V^ = A^-G, T 6j T} (22) 
V Vt = A t ~ -G tt 5j t < , 

ignoring the universal neutral current and gravitational potentials. The anti-neutrino 
potentials are the symmetric of these. 

In the following the numbers for the neutrino fluxes are extrapolated from the Su- 
pernovae neutrino luminosities obtained in with a 25M C model and in the generic 
model of [|nj. The energy fluxes follow this pattern: u^, v T , P^, v T have all the same lu- 



minosity say L Ut , and v e gets rapidly close to the v e luminosity after the initial v e burst 
(instant t — 0). From then to the supernova explosion (t ~ 0.45 sec), the average value of 
L Ve is around (4 — 8) x 10 52 ergs/sec and L Vt is about 2/3 to 3/4 of L Ue . After the explosion 
the spectrum is thermal and the luminosities L Vt ~ L Pe rs L Ve decrease along the time. 
At t = 1 sec, L Ue = (4 — 10) x 10 51 ergs/sec. The fluxes are obtained from the number of 
emitted particles per second F = j(4irr 2 ), and F = L/E by definition of average energy 
E. I take E = lOMeV for u e and 20MeV for v^v T . If the resonance is adiabatic i.e., 
all u e transform into v T and viceversa, the fluence transfer is 5F T = F Ue — F Vt . Then, in 
unities of 10 51 ergs/Mev/sec, bF^\ ~ 2 — 5 before the explosion and 5F^\ ~ 0.2 — 0.5 at 1 
sec time. In these unities, 

K M = -1.48G^ T ^xlO- 12 eV, (23) 
r io 

where G^ T = G^/Gp, Gf = ll,66TeV~ 2 is the Fermi constant and r w = r/10 10 cm. 

The standard charge current potential Vw is proportional to the electron density n e 
given by Y e ~ 1/2 times the nucleon density ~ p/m n . At large distances the mass density 
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p behaves as M/r 3 with a constant M = M 31 x 10 31 g, M 31 = 1 — 15 depending on the 



star \17\. Then 



V2G 



,. n e = 7.6Y e ^- x 1(T 13 eV . 
?io 



(24) 



Vjy decays necessarily faster than the goldston potentials which being proportional to the 
v fluxes inevitably overcome at large enough radii. That is a crucial point. The effective 
M 2 matrix [0 relevant for v e — oscillations is, up to a universal term, given by: 



( Ml = 2E{V i 



Vrj. 



2 e = Am 2 sin 29 



M 2 eil = Am 2 sin 29 



M 2 = Am 2 cos 29 



(25) 



An electron anti-neutrino is emitted with a negative potential Vp e = —Vw and a resonance 
would not be possible because Am 2 > 0, 9 < ir/4, by assumption. But, if the constant 
G^ T is negative, gets a negative goldston potential V p = G^ T 5j T and the resonance 
condition Ml = M 2 may be met. In the Fig. 1 it corresponds to the interception of the 
"e" and "/i" curves. 



10" 



T 1 — TTT 1 1 — TTT 1 1 — TTT 1 1 — TTT 1 1 — Til 




r (cm) 

FIG. 1. 2E(Vp e — Vp M ) is plotted in the solid and dotted curves with 5F§i = 1 and 3 respec- 

f 2 



tively. The dashed line stands for a particular value of M 2 = Am 2 cos 26 and the squares curve 



is the adiabatic boundary for M 2 tan 2 29 in resonance 2. 



The positive part of the 'potential' M 2 is plotted in Fig. 1 with values G^ T = — 1, 
5Fz,i = 1 or 3, Y e = 1/2, M31 = 4, and a neutrino energy of 20 MeV, the same for all curves. 
Take M 2 = 6 x 10~ 7 eV 2 for example. Its dashed line crosses the relative potential curve 
at two points: the first, r\, lies in the remarkably sharp transition from the —Vw < 
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dominated region to the goldston dominated one and very likely a non-adiabatic level 
crossing occurs there i.e., v e continues v e and continues v^. In the second point, r 2 , the 
potential is smooth enough (ex r -2 ) for adiabatic transitions. In general, the conversion 
probability is enhanced [|l2| , |l8f if M 2 tan 2 29 is larger than the rate f(r) = \E dV/Vdr\ 
at the resonance position. This function is plotted with squares in the figure for V oc r~ 2 . 
One can see that small mixing angle adiabatic conversions are possible for a wide range of 
neutrino masses. Using the resonance condition (neglecting V\y) one gets a lower bound 
on the mixing angle, 

M 2 tan^)>^^LxlO- 10 eV 2 , (26) 

while the maximum value of the v potential sets an upper limit to the square mass 
difference: 



Ml < 1.66Eio '^ T ~ xl(T 5 eV 2 , (27) 



where E 10 = E/lOMeV. 

The requirement of non-adiabatic crossing in the first sharp resonance puts an upper 
limit on M 2 tan 2 29. The slope of the potential on that point is equal (up to a numeric 
factor between 1 and 1/2) to the slope on the second resonance times the ratio (r2/Vi) 3 
and therefore that is the size of the window of interesting tan 2 29 values for a given M 2 . 
The net result is that for a wide range of parameters P e , i>a undergo one level crossing 



with maximum probability P c = 1 |T9| , pO| plus one adiabatic resonance. Their survival 
probabilities |i2|j20|| are therefore sin 2 9 and not cos 2 # (mixing with v T neglected). 

As far as u e , is concerned, the potentials and M 2 are symmetric to the anti- neutrino 
case. There is still one resonance position as in the standard model but now in the abrupt 
part of the potential. It results in a non-adiabatic crossing and the oscillation pattern is 
v e — > v T and v T — > u e , — > in contrast to the 'standard' one |IB|,I3|, u T — > v e — > 
Vfj, — > Ve, for the same mass mixing pattern. This does not change the observation 
predictions because v T have the same energy spectrum. Finally, the late v e 
oscillations alter the net flavour violation and modify the goldston fields out of the sphere 
where they take place but that does not create new resonances. 

This all picture may relate to solar neutrinos in two different ways: 1) v e — > v T in the 
Sun with 10~ 4 eV 2 > Am 2 T > Am^. Then, v e — > v T [y T — > is e ) triggers the appearance 
of goldston fields in the supernova (only out of the sphere where v e u T takes place 
in particular, if Am 2 T ~ 10 _5 eV 2 , r >2x 10 10 cm) and eventually v e <-> occurs. 2) 
Am 2 T lies above 10~ 4 eV 2 and it is v e — > that takes place in the Sun. Then, resonant 
v e *-* transitions occur in the supernova with the kinematic parameters proper of solar 
neutrinos. The solar non-adiabatic solution [|I9|-|2T| with parameters M 2 < eV 2 and 



M 2 tan 2 29 ~ 4 x eV 2 satisfies the desired condition if the neutrino fluxes are high 
enough (compare the dotted 5F 5 i = 3 curve and the solid one, 5F 5 i = 1). 

This points to a distinctive feature of that phenomenon: the strong dependence on the 
neutrino fluxes (third power in the bound of Eq. (P7|)) should manifest by the detection 
of P e s in the first moments of supernova observation with a energy spectrum proper of 
Pfj, (25 MeV instead of 15 MeV) that is not maintained later when the fluxes drop down. 
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A crucial condition is of course that the scale of Lepton numbers symmetry breaking lies 
at or below the Fermi scale and that the goldston potentials have the right sign. The 
observation of the described correlation between oscillations and flux magnitudes would 
provide a measurement of the breaking scale. 

The analyses of the SN1987A events \22\ have either disfavoured the occurrence of v e 
oscillations by comparing the best-fit energy spectrum with the current predictions for 



v e and v„ spectra P3J or put a maximum limit of 35% for the v e <-> z/„ permutation 



probability p4| . However, the oscillations driven by goldston fields, if they happen at all, 
only occur in the periods of largest fluxes. If they occur in the first 1/2 second or so of 
emission, the resultant hotter spectrum may manifest as a short artificial cooling time 
scale on the top of the real neutrino spectrum evolution. The fits to the SN1987A time 



evolution indeed exhibit two time scales, one of about 4.0 s and another < 1.0 s [53]. This 
seems consistent with an intrinsic two phase cooling process |25j nevertheless, it remains 
as a provocative piece of evidence. 

Another point that should be kept in mind is that v e oscillations appear |26j to 

be consistent or even favoured by the data when the analysis is made model independent. 
This may call for a spectrum of energies below those predicted so far but, as has been 
remarked recently P7| , this may be obtained if non-elastic neutrino scatterings in the 
transport sphere are taken in consideration. 

A typical signature of the assumed v mixing pattern is the non-observation of the 
initial v e burst. It turns more subtle however when one approaches the non-adiabatic 
frontier because the level crossing probability [|l^,^] depends exponentially on the inverse 
of the neutrino energy. Thus, almost-adiabatic transitions for 10 MeV neutrinos may be 
consistent with the survival of the highest energy neutrinos present in the electron capture 
v e burst. The survival of 32 MeV z/ e s was studied in [^| and in particular, its consistency 
with the solar non-adiabatic solution. 



IV. CONCLUSIONS 

The main result of this paper is that despite the fact that true Goldstone bosons only 
have derivative couplings with matter particles, one can prove that coherent Goldstone 
fields develop whenever the matter charges they are associated with are not conserved 
in macroscopic scales. The equations of motion of the Goldstone bosons are equivalent 
to the conservation law of the Noether currents associated to nonlinear symmetries of 
the type — > + a, these currents are a sum of a Goldstone boson piece and a matter 
current, and, whenever the quantum numbers of the matter particles are not conserved, 
the Goldstone fields develop compensating terms. Examples are where this may take 
place in a macroscopic proportion and as a result, the matter particles become subject to 
long range forces. 

In part because the potential energies are inversely proportional to the square sym- 
metry breaking scale and on the other hand, proportional to the time rates of particle 
reactions, the situations of interest are special, but not unique! Here I exercised a scenario 
where Lepton flavour violation is realized through neutrino oscillations in a Supernova 
with large enough intensity to alter the v potentials and oscillation patterns, provided 
that the symmetry breaking scale is less than 1 TeV. The caracteristic feature of those 
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kind of oscillations is their correlation with the periods of largest v fluxes thus providing 
a challenge to better understand the time evolution of the neutrino signal and the flavour 
dependent energy spectra as well. It seems worth to explore other physical systems and 
symmetries. 
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